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Abstract. In this paper we study the evolution of the wave function with the system 

size in a locally periodic structure. In particular we analyse the dependence of the wave 

^> I function with the number of unit cells, which also reflects information about its spatial 

r \ • behaviour in the system. We reduce the problem to a nonlinear map and find an 

• , equivalence of its energy regions of single periodicity and of weak chaos, with the 

' forbidden and allowed bands of the fully periodic system, respectively. At finite size 

the wave function decays exponentially with system size, as well as in space, when the 

energy lies inside a region of single periodicity, while for energies in the weak chaotic 

region never decays. At the transition between those regions the wave function still 

^ , decays but in a (/-exponential form; we found that the decay length is a half of the 

'^ ' mean free path, which is larger than the lattice constant. 
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1. Introduction 

Crystalline materials occupy a special place in the solid state physics. Despite that in 
real life the majority of the solids are non crystalline, crystals help to understand a lot 
of the properties of the solid matter. At the macroscopic scale crystalline solids are 
considered as infinite periodic systems, which are well described by the band theory [1] . 
However, the propagation of electrons through a locally periodic system, which consists 
of a finite number N of repeating elements [2], has been of great interest due to the 
practical applications in designing artificial materials with specific features. That is the 
case of layered periodic structures or finite superlattices [31 SJ [5] , and man-made devices 
fabricated with optical lattices to model condensed matter systems [6l [TJ [H [9l ITU] . 
Furthermore, the applicability of locally periodic structures cover a wide range that 
includes structurally chiral materials [H]; microwave [I2], photonic [iHl [HI [15], and 
phononic [HI |T7] crystals, as well as elastic or acoustic systems [HI |T9l [20] . 

The study of locally periodic structures is mainly concerned with the band 
formation for several cases of fixed N, where each one is analysed separately to look 
for the emergence of the allowed and forbidden bands [2ll [22], [23] [24]. In all cases, 
precursors of the band structure of the corresponding fully periodic system, which is 
almost formed when N is large enough, can be observed [2|. An interesting question 
that immediately arises is how or when the properties of a fully periodic systems are 
reached. Of course, the standard method of band theory for crystals is not applicable for 
finite systems. Therefore, to answer this question it is necessary to study the evolution 
of the system as A^ increases. This procedure allows to get additional information about 
the band formation; that is, the evolution to the infinite system should not be the same 
for energies in the gap than those in an allowed band of the fully periodic system. 

Along this line, the electronic transport on a double Cayley tree was studied as a 
function of the generation, which plays the same role as A^ [25]. There, a remarkable 
equivalence between the electronic transport and the dynamics of an intermittent low- 
dimensional nonlinear map has been exhibited; the map presents regions of weak chaos 
and of single periodicity, whose Lyapunov exponents are zero and negative, which 
indicates the conducting and insulating phases, respectively [26]. The conductance 
oscillates with the generation in the weakly chaotic attractors of the map, indicating 
conducting states. In the attractors of single periodicity, the conductance decays 
exponentially as is typical for insulating states, but at the transition the intermittency 
of the map makes the conductance to reach the insulating phase in a g-exponential form. 
Therefore, there exist typical length scales on the different energy regions that remain 
to be understood. 

A scaling analysis of the conductance is one way to obtain information about the 
metallic or insulating behaviour of a system [27]. But, the conductance is a global 
quantity and not much can be said about the local character in the system. Therefore, 
it is important to study the wave function directly, whose spatial behaviour should 
reflect whether a system behaves as a conductor or not. 
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The purpose of the present paper is to analyse the evolution of the wave function 
with A^ and once this is done, look at its spatial behaviour in a locally periodic structure 
in one dimension. This study allows to determine the nature of the several states in 
the system, to answer the question posted above, and to find a physical interpretation 
of the typical lengths scales in the problem. In order to do this, we consider a locally 
periodic potential that is formed by a sequence of individual scattering potentials. For 
our purpose, it is not necessary to consider an open system on both edges. Therefore, 
without any loss of generality, we avoid unnecessary complications and restrict the 
system to have only one entrance. We solve the problem analytically for arbitrary 
individual potentials using the scattering matrix formalism. Following the same ideas 
of reference [25j we reduce the problem to the analysis of the dynamics of an intermittent 
low-dimensional nonlinear map. 

We organize the paper as follows. In the next section we present the scattering 
approach to our linear chain of scatterers. In the same section, the recursive relation 
satisfied by the scattering matrix, when an individual scatterer is added, is reduced to 
a nonlinear map. The scaling analysis and the spatial dependence of the wave function 
is presented in section |3l We present our conclusions in section HI 

2. A serial structure of quantum scatterers 

2.1. Scattering approach 

The quantum system that we consider is a one- dimensional locally periodic whose 
separation between adjacent scatterers is a, that we will refer to it as the lattice constant. 
This chain consists of A^ identical scatterers, each one described by the potential Vb{x) 
of arbitrary shape and range b, as shown in figure [U The unit cell consists of a free 
region of size a — b plus the potential of range b, except in the first one where the size 
of the free region is a — b/2. Our chain lies in the semispace at the right of the origin, 
bounded on the left (x = 0) by a potential step of high Vq ^ E, where E is the energy 
of the quantum particle; the right side of the chain remains open. We assume that each 
scatterer is described by a unitary scattering matrix Sb that has the following structure: 

where vj, (r^) and t), (t|,) are the reflection and transmission amplitudes for incidence on 
the left (right) of the potential, respectively. 

As we will see below, the evolution of several properties of the system with 
the number of scatterers can be obtained from the scattering matrix of the system. 
Therefore, we write the scattering matrix S^, that describes the system with N 
scatterers, in terms of the one with A — 1 of them, Sn-i, by adding an individual 
scatterer using the combination rule of scattering matrices. The combination of Sn-i 
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Figure 1. A serial structure of identical scatterers, each one represented by a potential 
Vb of arbitrary shape and range &, whose associated scattering matrix is Sb- The chain 
is bounded on the left side by a potential step and opened on the right side. 



and Sh gives the following recursive relation, namely 

, 1 

< + h 
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-SN-lth 



b' 
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where k = \f2MEjW^ with M the mass of the particle. The initial condition for this 
recursive relation is the scattering matrix Sq associated to the scattering due to the 
potential step but measured at x = 6/2. 



2.2. Reduction to a nonlinear map 

Since S^^ relates the amplitude of the outgoing plane wave to the amplitude of the 
incoming one to the system, it is a 1 x 1 unitary matrix that can be parameterized just 
by a phase, 9iy, as 

^ = e^'". (3) 

I can be seen as a one-dimensional nonlinear map. That is. 



S, 



The recursive relation 
On = fidN-i), where 

f{ON-l) 



-6n-i + 2 arctan 



lm{r'f^al + abe'^^-^) 



Re (riaf + ate 



iOa 



(4) 



modulo TT. Here, ah = tbe*"'^/^e*''*^°~*\ e^'^' = t[/tb, and Re(a') and Im(a') denote the real 
and imaginary parts of a'. Because the reflected wave at the potential step acquires an 
additional phase 6'step, then Sq = e*^° where 9q = 9step + kb, which is the initial condition 
for the map. 

Although the nonlinear map (jl]) is valid for arbitrary individual potentials, it is 
necessary to consider a particular case to look for the characteristics of the map. Let us 
assume for a moment that an individual potential is just a delta potential of intensity v 
and null range, 6 = 0. In this case, rt, = r'j^ = —u/{u — 2ik) and t^ = t'^^ = —2ik/{u — 2ik), 
where u = 2Mv/h'^ [28]. The corresponding bifurcation diagram is shown in figure [2] (a), 
from which we observe ergodic windows between windows of periodicity one. The width 
of an ergodic window depends on the potential intensity, as well as on the energy; it 
becomes wider when ka increases (when referring to the delta potential we will speak 
of ka instead of k), as it is expected since for higher energies, smaller the effect of the 
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Figure 2. (a) Bifurcation diagram for a chain of delta potentials with ua — 10. We 
plot only the last thirty iterations of A^ = 1000 for the phase as a function of ka 
starting with an initial condition ^o — ^- We also plot the analytical results given 
by Eqs. (|12|) and (fT7|) but they are indistinguishable from the numerics, (b) Finite N 
Lyapunov exponent Ai(A^) as a function of ka for N — 1000. Theoretical result Ai 
given by (P^ is also plotted; it is indistinguishable from Ai(A^) for very large N. Inset: 
Ai for N = 20. 



potential. In the ergodic window 6^ fluctuates, while in the window of periodicity one 
it reaches a fixed value when A^ is very large. 

Figure |2] (a) suggests that, in a window of periodicity one, (E]) reaches a fixed point 
solution for S^ of the form shown in (|3]), for very large A^. If we look for it, we find 
a stable and an unstable fixed point solutions, e'^± [see figure |3] (a)]. In an ergodic 
window, (|2]) does not have solution of the form ([3]), but of the form w± = |w±|e*^, being 
9 the value around which 9n fluctuates with an invariant density. We will see below 
that these fixed point solutions are marginally stable [29]. The stable and marginally 
stable fixed point solutions can be summarized as 

3*^+ for kr' < k < kr^ 



S^ 



w± 



for 
for 



rCf "^ rv "^ rvf^f 



(5) 



kr' < k < k, 



Cm+l 



where k^^ and k^'^ denote the critical values of k on the left and right edges of the 
corresponding ergodic window [see flTO|l below]; 

1 



je± 






±^/ {Re abY — \t\\'^ + Hmaft 



for |ReQ;b(A;)| > |tfe(fc)P and 



w± 



I* 
n ah 



±v/|tfe|^-(Reafe)2 + Ima6 



(6) 



(7) 



for 



\Reab{k)\ < |4(A;)p 



(8) 
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such that 

tan^ = ^^^. (9) 

The equahty in (J8]) determines the critical values fee™ and fcc^; that is 

\ReabiK)\ = MK)\', (10) 

where kc denotes kc^ or kc'^. At a critical attractor, 6 takes the value 6{kc) = 6c, where 

_ lm[tr',{kc)al{kc)] 
'''''''- Re[tri{kc)al{kc)y ^^^> 

The condition (J8]) is equivalent to that for allowed bands in the infinite linear chain 
of scatterers [SU]. Therefore, we have found a correspondence between windows of single 
periodicity (or chaotic) and forbidden (or allowed) bands in the limit A^ — )■ oo, in a 
similar way as happens for the double Cayley tree studied in Ref. [25]. Equation (ITOl) 
defines the left and right edges of the chaotic windows. 

The dynamics of the map is expected to be different on each type of windows of the 
bifurcation diagram. To observe the trajectories on each region of k we need to back to 
the particular case of delta potentials. For this case, (|6]) gives 

2kn 

where 



x(ka) = \ { cos ka H sin ka 

^ ' y \ 2ka 



11(2 

y{ka) = sin/ca — cos ka (13) 

jdrCd 

for |cosA;a + {ua/2ka) sinA;a| > 1, while ([7]) gives 

2ka 

w± = €'''''' [±x'{ka)+y{ka)], (14) 

ua 

where 



x'{ka) = \jl — [ cos ka H — — sin ka ) (15) 



for 



ua 
cos ka H — ; — sm ka 
2ka 



< 1, (16) 

which is the condition for allowed bands in the Kronig- Penney model [31]. From ( 1T¥1) it 
is easy to see that the phase of w± is given by [see (|9])] 

= -ka + {m+l)Tc, m = 0, 1,... (17) 

The band edges are obtained from the equality in ( TT6l) . namely 

tan ^"^^ -[ W2A;c„a, form odd 
2 I ~2kc^a/ua, form even 

and kc'^a = mvr, for m odd and even. At the band edges, 0{kc^a) = 9c^, with 
6c^ = —kc^a + {m + 1)tt. For the first chaotic window {m = 1), k^^a = 2.28445. . . 
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Figure 3. Dynamics of tlie map (|4]) applied to tlie delta potential, where ua = 10 and 
ka taking the values (a) 2, (b) 2.2846, (c) 2.8, and (d) 3.14145. Each panel shows the 
trajectories for initial conditions 9q — 2 and 5. 



and k^ia = tt, 9^^ = 3.99873 . . ., and 6c' = vr. Equations (IT2l) and flTTl) coincide with the 
numerical results in both sides and inside of the first chaotic window, as it is shown in 
figure [2] (a). 

The trajectories of the map for this example are shown in figure |3] for some values 
of ka. In panel (a) ka = 2, ka is inside the window of period one; we observe that 
all trajectories (we show only two) converge to a fixed point as N ^ oo, which 
corresponds to 6+ ~ 3.61. Panels (b) and (d) show intermittent trajectories at the 
tangent bifurcations when ka = 2.2846 and 3.14145, very close to the transitions from 
the chaotic side. In panel (c), deep inside the chaotic window, the trajectories never 
converge. Therefore, we observe that the behaviour of the trajectories for two initial 
conditions strongly depends on the specific region to which ka belongs. The dynamics 
of the map can be analysed by means of the sensitivity to initial conditions, something 
that we do next. 



2.3. Sensitivity to initial conditions 

The dynamics of the nonlinear map 01]) is characterized by the sensitivity to initial 
conditions. For finite A^ it is defined by [25] 

d6j\f 



■^N 



dOn 



,AfAi(Af) 



(19) 



where 6q is the initial condition and Ai(A^) is the finite A^ Lyapunov exponent, whose 
dependence on A^ is shown as an argument (the subscript 1 will be clear below). For 
the map (|1]), E^ satisfies the following recursive relation: 

|4 



,NAi{N) 



U 



,(Af-l)Ai(Af-l) 



(20) 



N~l\ 
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In figure [2] (b) we plot Ai(iV) obtained from (1201) for the delta potentials as a 
function of ka, for A^ = 1000. We observe that Ai(A^) is negative in the windows of 
period one, indicating that S^r decays exponentially with A^ when A^ is very large. This 
means that any trajectory in the map converges rapidly to a fixed point, as happens 
in figure [3] (a). In the chaotic windows, Ai(iV) goes to zero in the limit N — )■ cx3. 
There, S^r does not depend on A^ and nothing can be said about the convergence of any 
trajectory; this situation corresponds to figure [3] (c). For a finite number of iterations 
Ai(A^) is still negative in the windows of period one, but it oscillates around zero in 
the chaotic windows. This can be seen in the inset of figure [21 the amplitude of those 
oscillations tend to zero as A^ increases; this is a signal of weak chaos [261 [32] . Here, we 
are interested in the behaviour of Sjv with A^. 

With this evidence, we can assume that, in the limit A^ — > oo, Ai(A^) — !■ Ai and 
Sat — >■ ^7v, where ^n is the sensitivity to initial conditions defined by 

eiv = e^^S (21) 

where Ai is the Lyapunov exponent of the map, which is given by 

A, = ln ^^ .. (22) 

Koil + abSool 

For one of the two roots expressed in ([6]), which are valid for fc in a window of period 
one, Ai is positive. This means that the fixed point solution is unstable. The second 
solution is a stable one since Ai is negative. It is the last one, the only that appears in 
figure [2] (a). In the chaotic windows there are two values of Ai that correspond to the two 
solutions w±, one positive and another negative, both at the same distance from zero. 
We find that the Lyapunov exponent that agrees with Ai(A^), which has been obtained 
iteratively for very large A^, is the average of those values, which is zero. In figure [2] (b) 
we compare Ai of f[22|) with Ai(A^) calculated iteratively by means of f[20l) . for the delta 
potentials. For A^ = 1000 we observe that both results are indistinguishable. Therefore, 
(pTj) says that C,n decays exponentially with A^ for k in the windows of period one and 
remains constant in the chaotic windows. These results are verified in panels (a) and 
(c) of figure [U where we compare (121]) with the numerical calculations for the delta 
potential. 

The behaviour of the sensitivity to initial conditions is very different at the 
critical attractors, since there an anomalous dynamics occurs due to the tangent 
bifurcation [33l[34]. A critical attractor is located at the point (fee, 6'c), where kc denotes 
kc^ or kc'^ and 9c, 9c^ or 9c'^. If we make an expansion of 9n close to 9c the result is 
given by 

9n~9c = {9n-i -9c)+u {9n-i - 9cy + ■ ■ ■ , (23) 

where z = 2 and u = =F|r[,(A;c)p. From known properties of this nonlinearity of the 
tangent bifurcation, the sensitivity obeys a g-exponential law for large A^ [3lj. That is, 
Sat^oo = ^N, where 

^^ oc e"!'^ = [1 - (g - l)A^A,]-^/(''-^) , (24) 
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Figure 4. The continuous lines represent the resuh for the sensitivity to initial 
conditions for finite iV, S^r, as a function of N for the chain of delta potentials with 
ua = 10 and ka equal to (a) 2, (b) 2.28445 . . ., (c) 2.8, and (d) 3.14145. The dashed 
lines are the theoretical results given by (|21[) for (a) and (c), and (|24|) for (b). In (d) 
^AT remains constant with N . 



with g = 1 — 1/z = 3/2 and A3/2 = zu = =F2|r^(/cc)p. The minus and plus signs 
correspond to trajectories at the left and right (right and left) of the point of tangency 
6'c = (^cm (^c = (^c'^)y respectively; that is, ^^ decays with A^ with a power law when 
On-1 — ^cm < i^N-i ~^c;„ > 0) and grows faster than exponential when 9n-i — ^c^ > 

In panels (b) and (d) of figure H] we show the behaviour of Sjy with A^ at the edges 
kc^a and kc'a, obtained from the numerical calculation for the delta potential, and 
compare them with (!24|) . In (b) it is clear that Sjv behaves as C,n of fl24|) at ka = k^j^a 
for very large N. That is not the case at the right edge ka = k^' a where S^r remains 
constant with N . This pathological behaviour is because fee' « corresponds to a resonance 
and the nodes the delta potentials become invisible. 



3. Behaviour of the wave function 



3.1. Evolution of the wave function with the system size 

The wave function in the region between the individual potentials can be written as a 
superposition of plane waves traveling to the left and right. If we normalize it, in such 
a way that the amplitude after the last scatterer is one, the square modulus of the wave 



function in the region between the scatterers n and n+1, labeled by n {n 
can be written as 



0,1, 






Wi 



XI 



NAgiN) 

fg 

nAq(n) 



cos^ [k{x -na- 6/2) + ft„/2] 



iV), 



(25) 



where 5' = 1 for k in the windows of single period and of weak chaos; q = 3/2 at the 
transition from the chaotic side. We see from 025p that at the position x, the amplitude 



Evolution with size in a locally periodic system: scattering and deterministic maps 10 

of the wave function depends on A^ only through the numerator. This means that we 
can replace A^ by A^ — 1 to find the wave function at the same position for a chain made 
of A^ — 1 scatterers. Hence, the following recursive relation is satisfied by the square 
modulus of the wave function: 

NAqiN) 

IV^f Ha:)!' = J\,„<._u \^L''-'\x)\\ (26) 



(7V-1)A,(JV-1) 



An iteration process lead us to 

|V^fna:)r = er'(^)|^i°)(x)f, (27) 

where ipn (x) is the wave function at x in the absence of any scatterer. 

The factor in front of the right hand side of (1271) is just the sensitivity to initial 
conditions for finite A^ [see f lT9|) ]. Therefore, the evolution of \-ipn (x)]"^ with A^ is given 
by Sjv, as shown in figure HI This evolution resemble the behaviour of the conductance 
in a double Cayley tree with system size [25]. For energies in a window of period 
one, \ijjn (a^)P decays exponentially with A^, for A^ :» 1, with a typical decay length 
which we identify with a localization length due to the similar scaling behaviour of the 
conductance with A^ [25]; this localization length is given by 

^1 = ih- (28) 

When I All > 1, on the one hand, the localization length is smaller than the lattice 
constant; this occurs far from the transition to the chaotic side; close to the transition 
Ci > a. On the other hand, for energies in the chaotic window, {ipn (a^)P never decays 
but oscillates as A^ increases; there, Ai = such that a localization length can not be 
defined. However, at the transition from the chaotic side the wave function shows a 
power law decay with A^ [see figure Hl^b)]. The typical decay length, that we identify 
with a localization length too, is given by 

What is very interesting here is that the term on the right hand side of the second 
equality is i/2, where 

i= ^, (30) 

which is a definition of the mean free path [35l 136] . This equation means that the mean 
free path is larger than the lattice constant. This result is in agreement with the one 
obtained from Esaki and Tsu, where i ~ 3a, in a supperlattice [37]. For our particular 
case of a delta potential with ua = 10, £ ~ 1.21a for k^^a = 2.28445 . . ., and i ~ 1.4a 
for kc' a = IT. The behaviour of the wave function in the space can let us to understand 
the nature of the states on the different energy regions. 
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Figure 5. Square of the wave function as a function of x for the chain of deha 
potentials with ua = 10. In (a) we compare p2|) (continuous hne) with an exponential 
with localization length given by ([SS)) (dashed line), for ka — 2. In (b) we compare 
(p3| (continuous line) with a g-exponential decay with localization length given by (|29|) 
(dashed line), for ka — 2.28445 . . .. (c) For ka = 2.8 the wave function is extended. In 
(c) the system is invisible at a resonance at ka — n. 



3.2. Spatial behaviour of the wave function 

In the limit of very large A^, Ag(A^) — )■ Ag and the squared modulus of the wave function 
given by ( 125|1 can be written as 



\t 



Wi 



X) 



nXn 



cos^ [k{x -na- 6/2) + 9^/2] . 



(31) 



For k far from the critical attractors, q = 1 and the ordinary exponential function is 
recovered. In this case the spatial behaviour of the wave function is given by 



hA. 



(N), 



X) 



-{N-n)a/a ^Qg2 ^^^^ _ ^^ _ ^^3) + ^„/2] . 



(32) 



This equation means that when the energy is in a window of period one, the wave 
function decreases exponentially in space from the boundary {n = N) to inside of the 
system (n < A^). In this case the quantum particle is spread over few Ci's. It is in this 
sense that we interpret Ci as a localization length. This localization length is smaller 
that the period of our locally periodic structure, a, except very close to the transition 
region. This situation is illustrated in figure [5] (a) for the chain of delta potentials. 

For energies in the weakly chaotic windows Ci — ^ oo and the wave function becomes 
extended through the system, as can be seen in figure [5] (c). 

Near the transition, at the weakly chaotic attractors, (!3T|) can be written as 

2 



|^fH^)| 



(^1^) cos^[M.-.-V2) + .„/2], 



(33) 



where we used that A3/2 = =F'^/C3/2- The plus (minus) sign in fl33|) corresponds to 
the left (right) chaotic attractor. This equation clearly shows a power law decay from 
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the boundary. That is, the wave function is still localized but the localization is not 
exponential, but g-exponential. It is interesting to note that the localization length is 
C3/2 = ^/2. The factor 1/2 is due to the fact that the localization length is measured from 
the maximum of the wave function, which in this case is at the boundary of the system, 
while the mean free path implicitly assumes the width of a wave packet; according to 
Esaki and Tsu |37], the mean free path in a supperlattice is the uncertainty in the 
position. In figure [5] (b) we can observe the power law decay of the squared modulus of 
the wave function for the chain of delta potentials. Figure [5] (d) shows the behaviour 
of the squared modulus of the wave function at the attractor on the right side of the 
chaotic window. In this case, the system is invisible since k^' a = tt corresponds to a 
resonance and the wave function has nodes just at the delta potentials. 

4. Conclusions 

We considered a locally periodic structure in one dimension, which consists of a chain of 
potentials of arbitrary shape. We studied the evolution of the wave function of the chain 
when the system size increases. This procedure helped us to observe the behaviour of 
the wave function in space when the size of the system remains fixed. 

Since the system is not fully periodic, we can not use the traditional band theory. 
Instead of that we take advantage of a recursion relation of the scattering matrix, in 
terms of the number of scatterers, to reduce the problem to a nonlinear map. Through 
this method all information about the behaviour of the system was obtained from the 
dynamics of this map. In this way we could understand how and when the band theory 
of the fully periodic chain is reached by the knowledge of the type of evolution of the 
wave function with the system size. We found an equivalence between the periodic and 
weakly chaotic regions of the map and the forbidden and allowed bands, respectively. 
This equivalence is similar to that remarked in the literature for the conductance of a 
double Cayley tree, where it is null in the windows of period one, and oscillates in regions 
of weak chaos. Also, we found that the wave function at a given position scales with 
the system size in a similar way as the conductance does. In a window of period one, 
far from the transition to the chaotic window, the wave function decays exponentially 
with a typical length scale, which is smaller than the lattice constant. Very close to the 
transition from the window of period one, this typical length becomes larger than the 
lattice constant. This behaviour driven us to interpret the typical scale as a localization 
length. We corroborated this interpretation by the exponential localization of the wave 
function close to the boundary of the finite system. In the chaotic region the wave 
function does not decay and a localization length can not be defined. 

At the transition between periodic and weakly chaotic regions, but from the chaotic 
side, the wave function scales as a power law. This implies that it is still localized but 
with a g-exponential form. The localization length were found to be one half of the 
mean free path. We demonstrate analytically, that the mean free path is larger than 
the lattice constant, in complete agreement with a result found in the literature for a 
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supperlattice. 

Finally, it is worth mentioning that our development considers independent 
quantum particles such that classical waves can also be used. In particular, one- 
dimensional elastic systems are strong candidates to simulate our quantum system by 
means of elastic rods with narrow notches |19j . 
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